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Abstract 

The r decay rate into hadrons of invariant mass smaller than ^/sq 3> 
Aqcd can be calculated in QCD assuming global quark-hadron duality. It 
is shown that this assumption holds for so > 0.7 GeV 2 . From measure- 
ments of the hadronic mass distribution, the running coupling constant 
oe s (so) is extracted in the range 0.7 GeV 2 < s < m 2 T . At so = rn^, the 
result is a s (m%) = 0.329 ± 0.030. The running of a s is in good agreement 
with the QCD prediction. 



(Phys. Rev. Lett. 76 (1996) 3061) 



Ihe scale dependence ot coupling constants is one ot the key teatures ot renor- 
malizable quantum field theories. In QCD, the effective coupling constant a s (Q 2 ) 
is predicted to decrease with the momentum transfer Q 2 , a property referred to 
as asymptotic freedom This prediction has been tested by comparing data 
obtained from experiments operating at different energies ||; it has also been 
studied in single high-energy experiments at ep and pp colliders, where a large 
range in Q 2 can be probed simultaneously 0. Here we propose a test of the 
scale dependence of a s {Q 2 ) in the low-energy region 0.7 GeV 2 < Q 2 < m 2 . Our 
method is based on integrals of the invariant mass distribution in hadronic r 
decays. It provides a unique opportunity to test one of the most important pre- 
dictions of QCD in a single experiment and at low energies, where the effect of 
the running of a s is most pronounced. 

We shall consider the r decay rate into hadrons of invariant mass squared 
smaller than sq, normalized to the leptonic decay rate: 

So 

r(r -> v r + hadrons; s ha d < s ) f dR T (s) 

Rr (s ) = = / ds — , 1 

F(t — > u T e v e ) J as 

where dR T /ds is the inclusive hadronic spectrum. As long as so ^> Aq CD , the 
quantity R t (sq) can be calculated in QCD using the Operator Product Expansion 
(OPE) 0, |5|. Applying the OPE in the physical region assumes global quark- 
hadron duality, i.e. that decay rates admit a QCD description after a "smearing" 
over a sufficiently wide energy interval has been performed ||, which in the 
present case is provided by the integration over the range < s < s . The 
question of how accurate this assumption is and for what values of Sq it applies 
is a phenomenological one; it cannot be answered yet from theoretical grounds. 
Below, we shall investigate this question, comparing data with theoretical pre- 
dictions based on the duality assumption. A similar test of duality has been 
performed in Ref. J7|, using data on the e + e~ — ► hadrons cross section. 

The t decay rate into hadrons can be written in terms of moments M.^ °f the 
absorptive part of current-current correlation functions of angular momentum J 
|§ §]. The quantity R t (sq) is given by 



R T (s ) = ^Mt\s ) - 2 (4) M?(s ) + (4) M£\so) 
mi \miJ \mij 
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+ ^M 0) (so)-2(4) Mf\s a ) + 2 -(^] Mf\s ), (2) 



mi \mij o \m 



where Sew — 1.0194 accounts for electroweak radiative corrections [ID|. The 
moments can be written as contour integrals along a circle of radius sq in the 
complex plane. Since the only large mass scale in these integrals is Sq, the OPE 
provides an expansion in powers of 1/sq: 

M J) (*o) =-Ml 1) [«,( S o)U^=i + £^ J) [« s (3o)]%^. (3) 

71=1 S 
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Ihe leading term is given by perturbation theory alone, terms suppressed by 
powers of 1/sq consist of perturbative coefficients multiplying dimensionful 
parameters (0 2n ), such as quark masses or vacuum condensates ||. This is 
how nonperturbative effects are incorporated. There is no leading term for the 
moments with J = 0, which vanish in the chiral limit and are thus proportional to 
powers of the light quark masses. For the moments with J = 1, the perturbative 
contribution is 

((1)U l, M _ 1 , ^ Mk) fUs(So) 



a, 



(so)) 



pert 



n=l 



7T 



where a s (so) is defined in the MS renormalization scheme, 
three coefficients are given by M, |[ 
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1, and the next 
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fc + 1 (k + 1) 2 (k + 1) 



The coefficient Ka appears in the perturbative expansion of the Adler function 
and is not known exactly. An estimate using the principle of minimal sensitivity 
and the effective charge approach |Tl| gives Ka ~ 27.5 |12|]. We shall use this 
result in our analysis. The error due to the truncation of the perturbation series 
in (Q) is of the order of the last term included. It can also be estimated by 
summing a subset of corrections to all orders in perturbation theory. Such a 
class of corrections is provided by the renormalon chains fl3[l , which are the 
terms of order /?Q _1 a", where /?o is the first coefficient of the /5-function. For 
the case of the moments, the resummation of these terms has been discussed 
in Refs. 
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Below, we shall take fixed-order perturbation theory as the 
nominal scheme and use the resummation of renormalon chains to estimate the 
perturbative uncertainty. 

The nonperturbative corrections in the OPE are proportional to the light 
quark masses or to vacuum condensates @]. We quote the power corrections for 
the sum of the moments contributing to R t (sq) in (0). The terms relevant to the 
numerical analysis are 
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where m s [s ) is the running strange-quark mass, \m q ip q ip q ) are the quark con- 
densates, and pa s {ijjip) 2 denotes the four-quark condensate. More detailed ex- 
pressions, which are used in our analysis, can be found elsewhere [|5], ||, §. At 
tree level, the powers of 1/sq appearing in the OPE of the moments in ([3]) con- 
spire with the powers of So/m 2 ., which multiply the moments in @, so that the 
nonperturbative corrections to R t (sq) are suppressed by powers of 1/m 2 .. This 
is no longer the case if radiative corrections to the coefficients of the vacuum 
condensates are taken into account, but the corresponding effects are very small. 
As a consequence, the power corrections to R t (sq) remain small down to rather 
low values of s ; using standard values of the QCD parameters (which we take 
from Ref. 0) we find —(1.4 ± 0.5)% for the right-hand side of @ at So = rn%., 
and —(1.5 ± 0.5)% at so = 1 GeV 2 . This observation, together with the fact 
that the perturbative contributions are known to high order, guarantees a good 
convergence of the OPE down to low energy scales. 

To extract the quantity R t (sq), we use the spectra of the hadronic mass dis- 
tribution reported by the CLEO and ALEPH Collaborations [15|, [U| (see Fig. [I]). 
To obtain dR T /ds, we multiply the normalized distributions by the world average 
R T = R T {m 2 .) = 3.642±0.010 [|17]]. Not shown is the contribution from r — > hTv T 
with h~ = ir~ or K~ , which has a branching ratio of (11.77 ± 0.14)% We 
integrate these spectra over s, combine the results weighted by their statistical 
errors, and add the systematic errors, which we estimate by taking the differ- 
ence between the CLEO and ALEPH data.[] This is justified, since the dominant 
sources of systematic errors are different in the two analyses. The result is shown 
in Fig. [I]. It is represented as a band, since the errors in the R T (so) values are 
strongly correlated. The two curves show theoretical calculations of R t {sq) based 
on the OPE approach outlined above. The solid line is obtained using fixed-order 
perturbation theory to order af. The dashed line is obtained by adding to this 
a resummation of renormalon chains of order a 5 s and higher, using the results of 
Ref. 0. The value of a s {m 2 ) has been adjusted so as to fit the data at s = m 2 .. 
The central values obtained in the two schemes are a s (m 2 ) = 0.329 (fixed-order) 
and a s (m 2 ) = 0.309 (resummed). Their difference provides an estimate of the 
uncertainty due to unknown higher-order corrections, which is more conservative 
than that obtained by omitting the term of order a* in the fixed-order calcula- 
tion. Varying the values of the nonperturbative parameters within conservative 
limits changes a s (m 2 ) by up to 2%. Adding linearly the perturbative uncertainty 
(±0.020), the nonperturbative uncertainty (±0.006), and the experimental un- 
certainty (±0.004), we find 

a s (m 2 T ) = 0.329 ±0.030, 

a s (m 2 z ) = 0.119 ±0.004. (7) 

*As the ALEPH data are preliminary, this estimate may be taken with caution. However, 
since inclusive quantities such as R t (sq) do only probe gross features of the hadronic mass 
distribution, systematic errors play a minor role in our analysis. 
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.bor the sake ot completeness, we have translated our result into a value ot a s at 
the mass of the Z boson. 

The assumption of global quark-hadron duality can be tested by comparing 
the data for the quantity R T (s ) at values s < m 2 with the theoretical pre- 
diction ||. Given a s (m 2 ), the value of a s (s ) follows from the solution of the 
renormalizat ion-group equation 

/3K)=A^ + /3 1 (^) 2 + &(^) 3 + ..., (8) 

where /3o = 9, /3i = 64 and /?2 = 3863/6 are the first three coefficients of the 
/3-function, evaluated for rif = 3 light quark flavours. (The value of @2 is specific 
to the Ms scheme.) Whereas the theoretical uncertainties are a limiting factor in 
the determination of a s (m 2 ), they have little influence on the Sq dependence of 
R T (so). For the perturbative part of the calculation this is apparent from the good 
agreement of the two curves in Fig. [I], which refer to values of a s (m 2 ) that differ 
by 9%. Varying the values of the nonperturbative parameters has a negligible 
effect (~ 0.5% at Sq — 1 GeV 2 ) on the value of R T (so) for so < m 2 , since the only 
dependence on sq in (||) comes from the quark-mass corrections, which are known 
with higher accuracy than the vacuum condensates. Hence, the Sq dependence 
of R t (sq) is predicted essentially without free parameters, and the comparison 
between theory and experiment provides a direct test of quark-hadron duality. 
We find good agreement over the entire range 0.7 GeV 2 < so < m 2 , indicating 
that in r decays duality holds as soon as the integral over the hadronic mass 
distribution includes the p resonance peak. The onset of duality happens almost 
instantaneously, in accordance with general expectations. 

From now on, we shall rely on this behaviour and assume that for s > 
0.7 GeV 2 possible violations of duality can be neglected. We then turn to the 
main focus of our study: a test, at low energies, of the QCD prediction (||) 
for the running of the coupling constant, r decays are an ideal place to study 
this phenomenon since the value of a s (so) changes by a factor 2 over the region 
where duality holds, which is about the same change as in the region between 5 
and 100 GeV. From the quantity R T (so) shown in Fig. [I], we extract a s (so) as 
a function of so by fitting to the data the theoretical prediction obtained from 
(|D, (||)-(|)- The result, including experimental errors only, is represented by the 
dark band in Fig. ^. Theoretical uncertainties arise from the truncation of the 
perturbation series and from the uncertainty in the values of the nonperturbative 
parameters. As discussed above, they affect the overall scale of the a s values 
(by about 8-10%), but have very little effect on the evolution of the coupling 
constant. The sum of the experimental and theoretical errors is represented by 
the light band. The dashed curve shows the QCD predictions for a s (s ) obtained 
at three- loop order, normalized to the central value of the data at s = m 2 . The 



observed scale dependence ol the running coupling constant is m good agreement 
with the QCD prediction. The small oscillation of the experimental band around 
the theoretical curve, which could be due to some deviations from duality in the 
ai region, are not significant given the precision of the data. 

To quantify this agreement, we extract from the data the /^-function that 
describes according to (|8|) the running of a s (so). Defining x = a s (so)/47r, we 
have 

A) + /3iX + (3 2 x 2 + ... . (9) 



a 2 (so) din so x 

We approximate the derivative da^/dlnso by a ratio of differences, Aa s /Alnso, 
for a set of so values chosen such that the differences Aa s are large enough to 
be significant given the errors in the measurement. For a s (so) in (|9]) we take the 
central value of each interval. We use the following s values: 0.75, 0.95, 1.35, 
2.06, and 3.16 GeV 2 , corresponding to four intervals of increasing width Ams , 
but constant Aa s ~ 0.075. The results are shown in Fig. ^. The estimate of the 
errors includes the theoretical uncertainties, the error due to the choice of finite 
intervals in a s , and the experimental errors, which in this case are the dominant 
ones. The curves in Fig. |3] show the QCD /^-function at one-, two- and three- loop 
order in perturbation theory. The data provide clear evidence for the running 
of the coupling constant. Moreover, they prefer a running that is stronger than 
predicted at one-loop order. Between the three curves, the one that shows the 
three-loop prediction provides the best description of the data. Performing a fit 
with the three-loop /3-function, where /3 = 9 and fti = 64 are kept fixed but the 
three-loop coefficient j3 2 is treated as a parameter, we find P^/Pf 1 = 1-6 ± 0.7. 
We believe that such an experimental determination of the /^-function beyond 
the leading order can at present be done only in r decays. (A high-precision 
measurement of R e + e - (s) in the region below the charmonium resonances would 
provide an alternative place for such a study.) At higher energies, the value of 
a s is too small to distinguish between the three curves in Fig. || measurements 
in the region Q ~ 100 GeV, for instance, correspond to values x ~ 0.01. 

In summary, we have presented a method to measure the running coupling 
constant a s (s) in the low-energy region 0.7 GeV 2 < s < m 2 using r decay data 
obtained in a single experiment. It provides a test of one of the key features of 
QCD in a region where the effect of the running of a s is most pronounced. The 
theoretical analysis is based on the OPE and the assumption of global quark- 
hadron duality. We have tested this assumption and find that it seems to hold if 
the t decay rate is integrated over an energy interval large enough to include the 
p resonance peak. Our analysis provides a test of QCD at scales comparable with 
the lowest ones attainable before (Q 2 ~ 2.5 GeV 2 in deep- inelastic scattering), 
and with higher precision than all other single measurements of the running to 
date. We have extracted for the first time the /3-function from data and find that 
it is in good agreement with the three- loop prediction of QCD. 
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Figure 1: The quantity R t (sq) extracted from the data on the hadronic 
mass distribution dR T /ds reported by the CLEO and ALEPH Collaborations 
15 , (small figure). The experimental result is represented as a band. The 



curves show the theoretical predictions (see text). 
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Figure 2: Values of a s (so) extracted from the data on R t (sq). The dark band 
represents the experimental errors, the light one the sum of the experimental 
and theoretical errors. The errors are strongly correlated. The dashed line 
shows the three-loop QCD prediction for the running coupling constant. 
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Figure 3: Experimental determination of the /3-function. The curves show the 
QCD prediction at one- loop (dash-dotted), two- loop (dashed) and three- loop 
(solid) order. 



10 



